EPJ manuscript No. 

(will be inserted by the editor) 



PITT PACC 12002, DESY 12-020, FR-PHENO-2010-030, 
MADPH-10-1562, EDINBURGH-2010-25 



IPHC-PHENO-10-03, IPPP/10/80, DCPT/10/160, 



Introducing an interface between FeynRules and WH IZARD 



Neil D. Christensen a ' 1 , Claude Duhr b ' 2 , Benjamin Fuks 0,3 , Jiirgen Renter' 



d,4,5,6 



and Christian Speckner' 



e,5 



o- 

(N; 
Or 

< 

(N 



(N 
> 

in 

(N 
m 

O 
o 



1 piTTsburgh Particle physics, Astrophysics and Cosmology Center, Department of Physics and Astronomy, University of 
Pittsburgh, Pittsburgh, PA 15260, USA 

2 Institute for Particle Physics Phenomenology, University of Durham, Durham, DH1 3LE, U.K. 

3 Institut Pluridisciplinaire Hubert Curien/Departement Recherches Subatomiques, Universite de Strasbourg/CNRS-IN2P3, 
23 Rue du Loess, F-67037 Strasbourg, France 

4 DESY Theory Group, Notkestrafie 85, D-22603 Hamburg, Germany 

5 Albert-Ludwigs-Universitat Freiburg, Physikalisches Institut, Hermann-Herder-Strafie 3, 79104 Freiburg, Germany 

6 University of Edinburgh, School of Physics and Astronomy, JCMB, The King's Buildings, Mayfield Road, Edinburgh EH9 
3JZ, U.K. 

February 2012 

Abstract. While Monte Carlo event generators like Whizard have become indispensable tools in studying 
the impact of new physics on collider observables over the last decades, the implementation of new models in 
such packages has remained a rather awkward and error-prone process. Recently, the FeynRules package 
was introduced which greatly simplifies this process by providing a single unified model format from which 
model implementations for many different Monte Carlo codes can be derived automatically. In this note, 
we present an interface which extends FeynRules to provide this functionality also for the Whizard 
package, thus making Whizard's strengths and performance easily available to model builders. 

PACS. 11.80.Gw Multichannel scattering 12.60.CnExtensions of electroweak gauge sector 12. 60. Fr Exten- 
sions of electroweak Higgs sector 



1 Introduction 



■ • ■ The Standard Model (SM) of particle physics provides 
. £h ! a successful description of all experimental high-energy 
■ data to date. However, despite its success, many funda- 
5^ ' mental questions remain unanswered, such as the origin 
5^ . of electroweak symmetry breaking, the nature of neutrino 
masses, the large hierarchy between the electroweak and 
the Planck scales and the origins of dark matter and the 
cosmological constant. Attempts to address these ques- 
tions have led to a wide range of new physics theories, 
most of them predicting new phenomena at the TeV scale. 
The Large Hadron Collider (LHC) at CERN is currently 
probing this scale and will hopefully make it possible to 
discover, constrain and/or exclude some of the theories 
beyond the Standard Model (BSM). 
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Monte Carlo simulations of the particle collisions to be 
observed at the LHC will play a key role in the exploration 
of the weak scale, both from the theoretical and experi- 
mental sides. While proper modelling of the strong inter- 
actions, including parton showering, fragmentation and 
hadronization, is essential in order to achieve a realistic 
description of the event distributions at the LHC, BSM 
signals are expected to show up predominantly in the un- 
derlying hard interaction. Therefore, a lot of effort has 
been put into the development of multi-purpose matrix- 
element and parton-level event generators such as Com- 
pHep/CalcHep QII3], MadGraph/MadEvent [1G3 
EMZMS], Sherpa [§] and Whizard [lOUllj . which allow, at 
least in principle, the generation of parton-level events for 
a large class of Lagrangian-based quantum field theory 
models. 

In practice, however, the implementation of a BSM 
model into any of these tools can be a tedious and error- 
prone task. Not only might a model contain hundreds, if 
not thousands, of interaction vertices which need to be 
encoded in a format suitable for the generator in ques- 
tion, but in addition, each code follows its own format 
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conventions which need to be respected. To improve this 
situation, a framework based on the FeynRules pack- 
age [12] , addressing not only the implementation but also 
the validation of new physics models in multi-purpose 
matrix-element and event generators has recently been 
proposed 43 and its virtue illustrated in the context of 
the programs CompHep/CalcHep, MadGraph/Mad- 
Event and Sherpa. 

In this paper we present a new interface from Feyn- 
Rules to the event generator Whizard. It follows the 
same approach as the already existing interfaces, and al- 
lows to export any model implemented in FeynRules 
directly to Whizard, extending thus the aforementioned 
framework to the Whizard package. The paper is orga- 
nized as follows: after giving a short overview of the Feyn- 
Rules and Whizard packages in Section [2l we discuss 
the usage, features and limitations of the new interface 
in Section [3] Finally, we give a few examples of using the 
interface in Section 2] before concluding in Section [SJ The 
Appendix contains an exhaustive list of the interface op- 
tions (Appendix |X| and a selection of numerical results 
from the interface validation (Appendix [Bjl . 



2 The Whizard and FeynRules packages 
2.1 FeynRules 

FeynRules is a Mathematical 5 ] package that allows to 
derive Feynman rules from any perturbative quantum field 
theory-based Lagrangian in an automated way. The in- 
put provided by the user is threefold and consists of the 
Lagrangian defining the model, together with the defi- 
nitions of all the particles and parameters that appear 
in the model. Up to now, the public release of Feyn- 
Rules supports scalar, vector, fermion (Dirac, Majorana 
and Weyl [14]) and spin- two fields, as well as Faddeev- 
Popov ghosts, while recently, superfields have also been 
included [TS]. Once this information is provided, Feyn- 
Rules can perform basic checks on the sanity of the im- 
plementation (hermiticity, normalization of the quadratic 
terms, . . . ), and finally computes all the interaction ver- 
tices associated with the model and store them in an in- 
ternal format for later processing. 

After the Feynman rules have been obtained, Feyn- 
Rules can export the interaction vertices to various matrix- 
element generators by means of interfaces provided by the 
package |13j . The current public release contains inter- 
faces to CalcHep/CompHep, FeynArts/FormCalc, 
GoSam [TB1IT7] . MadGraph/MadEvent (both 4 and 
5), Sherpa and the newly developed interface to Whizard 
which we present in this note. 

Since the matrix-element generators very often have 
color and/or Lorentz structures hardcoded, the different 
interfaces check whether all the vertices are compliant 
with the structures supported by the corresponding matrix- 
element generators, and discard them in the case they are 

1 Mathematica is a registered trademark of Wolfram Re- 
search Inc. 



not supported. The output of an interface consists of a 
set of files organized in a single directory which can be 
installed into the relevant matrix-element generator and 
used as any other built-in models. 

2.2 Whizard 

Whizard [11] is a program package for the automatic and 
efficient generation of unweighted parton-level events for 
multileg tree-level processes in both the Standard Model 
and a large number of BSM models. While the original ver- 
sion of Whizard (l.xx) was geared towards ILC physics 
by offering elaborate options for controlling the beam setup, 
the new version Whizard 2.x.x has been redesigned as 
a hadron collider tool vastly extending but also including 
the (ILC) features of version 1 (it now also includes both a 
fey-ordered and an analytic parton shower |18]L Although 
many QCD/NLO and SM background physics topics (c/. 
e.g. [20 , 14 ( 21 ( 19]) have been tackled using the package, 
its main focus lies on BSM physics and its collider phe- 
nomenology (cf. e.g. [55]). Specifically models focusing on 
alternative scenarios have been incorporated in Whizard 
and independently been tested and validated 22 , 24, 27 . 
28,29,30,31;, where specifically also new physics effects 
from higher-dimensional operators have been investigated. 
The new interface to FeynRules described in this paper 
makes such studies considerably easier. 

The program is not designed as a monolithic block, but 
instead shows a high degree of modularity including the, 
in principle, independent and stand-alone usable matrix- 
element generator O'Mega and the adaptive multi-channel 
Monte Carlo integrator Vamp as well as several dedicated 
sub-packages for lcpton collider physics, all of which are 
contained in an easy-to-use package. This design is shown 
schematically in Fig. [1] Whizard has been designed such 
that this structure is completely transparent to the user. 

For providing additional functionality like convolution 
with parton distributions, hadronization or support for 
additional event formats, Whizard can be linked against 
external packages like Lhapdf [35], Pythia [37] or also 
HepMC [38]. 



O'Mega: O'Mega [TU] is a generator for tree-level ma- 
trix elements which generates a symbolic representation 
of the scattering amplitude and translates it into For- 
tran 90 code. As this is a computer-algebraic process 
which is rather similar to the workflow of an optimizing 
compiler, the OCaml language has been chosen for the 
implementation of O'Mega. 

The algorithm is based on the algebraic paradigm of 
so-called Directed Acyclical Graphs (DAGs) which allow 
a re-usage of all already computed components. Inspired 
by the ideas of Helas [3j5] , amplitudes are constructed by 
building a representation of the matrix element as a tree 
of one-particle irreducible wave functions (lPOWs, Greens 
functions with all but one leg amputated) which is then 
translated into highly optimized Fortran 90 code to 
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WHIZARD core 
user interface, steering, phase space 



O'Mega 
matrix elements 



VAMP 
Monte-Carlo integration 



USER 

process setup, cuts, analysis definitions, etc. 

Fig. 1. General structure of the Whizard package. 



evaluate the matrix element by numerically fusing lPOWs. 
This approach avoids the repeated and redundant evalu- 
ation of subdiagrams present in any Feynman diagram 
based approach and brings the factorial growth in com- 
plexity of the latter down to an exponential behavior. In 
addition, the optimized structure of the amplitudes is es- 
pecially well-suited to cope with big cancellations present 
in theories with gauge invariances. 

The algebraic O'Mega algorithm bears similarities to 
the numerical recursive solution of Schwinger-Dyson equa- 
tions of motion implemented in the code Alpha [40] . In 
this program, however, the evaluation proceeds fully nu- 
merically and no symbolic representation of the amplitude 
is constructed. 

While there is no design limitation on the spin of the 
fields or on their interactions, O'Mega currently supports 
spin 0, ^, 1, | and 2 fields, most dimension three and four 
couplings between those (in particular all dimension four 
gauge interactions) and a small set of higher dimension 
operator^. 

Color treatment in O'Mega leverages the color flow 
[1T1H2"] decomposition and currently supports fields trans- 
forming either as singlets, triplets, antitriplets or octets. 
Only the vertex factor belonging to a predefined flow is 
implemented in the model definition, and the other flows 
are derived by O'Mega from the SU(3)c representations 
of the fields. 

Models are represented in O'Mega as OCaml mod- 
ules which contain symbolic definitions of the fields and 
coupling constants, vertex lists, translation functions which 
map the fields onto their Lorentz / SU(3)c representa- 
tions and a number of auxiliary maps which provide tex- 
tual representations of the fields and constants for inter- 
facing and code generation purposes. These modules are 
compiled and linked with the O'Mega framework to ob- 
tain executables which can be called to transform a pro- 
cess definition into Fortran code. 



2 A future revision of Whizard /O 'Mega will allow specify- 
ing the Lorentz structure of the couplings directly in the model 
file, thus alleviating the limitation to predefined structures and 
permitting a much larger set of couplings. 



Vamp: Vamp [13] is a multichannel extension of the Ve- 
gas [Hj algorithm. For all possible singularities in the 
integrand, suitable maps and integration channels are cho- 
sen which are then weighted and superimposed to build 
the phase space parameterization. Both grids and weights 
are modified in the adaption phase of the integration. 

The multichannel integration algorithm is implemented 
as a FORTRAN 90 library with the task of mapping out 
the integrand and finding suitable parameterizations being 
completely delegated to the calling program (Whizard 
core in this case). This makes the actual Vamp library 
completely agnostic of the model under consideration. 



Whizard core: With matrix element generation and Monte 
Carlo integration being outsourced into dedicated building 
blocks, it is the job of the program core to roll O'Mega 
and Vamp transparently into a package which can be op- 
erated without dealing with the components directly. 

To this end, the core provides a command-line based 
user interface which is used to specify the process, event 
generation and analysis setup and to control the runs of 
the event generator. In addition, it is responsible for all 
physics pieces not covered by the other two components, 
most importantly phase space generation and beam (par- 
ton densities) setup. 

The major task of the WHIZARD core is to generate 
the phase space and map out the integrand, for which 
it needs information on the model under consideration. 
In the legacy 1.x branch of Whizard, the core is imple- 
mented as a set of Fortran 90 modules, Perl glue for 
managing the code generation and Makefiles for steering 
the build process. New models have to be added to the 
build framework in order to be available to the user. 

The recently published 2.x branch of Whizard super- 
sedes the old versions and features a complete rewriting of 
the core as a single Fortran 2003 program, discarding 
the old Perl code. Models are loaded dynamically and 
searched for in a configurable search path, allowing the 
user to add new models without modifying the Whizard 
package itself. Moreover, the dedicated scripting language 
Sindarin has been introduced which steers all aspects of 
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the simulation and which replaces the set of input hies 
used in previous versions. 



3 Implementing new models in Whizard using 
FeynRules 

3.1 General strategy 

In this Section, we describe how to implement generic 
BSM models into Whizard in an efficient way by means of 
the newly developed interface from FeynRules to Whiz- 
ard. 

This new interface brings the implementation of new 
models into Whizard in line with the approach intro- 
duced in Ref. [T3], where its power to develop, implement 
and validate models has been demonstrated in the con- 
text of CalcHep/CompHep, MadGraph/MadEvent 
and Sherpa. This approach is built around FeynRules, 
a Mathematica package where any perturbative and lo- 
cal quantum held theory Lagrangian can be implemented 
and its Feynman rules obtained in an automated way. The 
interaction vertices can then be passed, together with all 
relevant information, to various matrix element generators 
through dedicated interfaces. 

The new FeynRules interface to Whizard allows 
the user to integrate any model implemented in Feyn- 
Rules into the Whizard/ O' Mega framework, although 
some technical restrictions apply and will be explained 
in the next section. The interface produces a set of hies 
that contain all the particle and parameter definitions of 
the model together with the interaction vertices. Together 
with the model hies, a framework is created which allows 
to communicate the new models to Whizard in a well 
defined way, after which step the model can be used ex- 
actly like the built-in ones. This specifically means that 
the user is not required to manually modify the code of 
Whizard/O'Mega, the models created by the interface 
can be used directly without any further user intervention. 



3.2 Features and restrictions of the Whizard - 
FeynRules interface 

In this Section we describe in detail the features and re- 
strictions of the new interface from FeynRules to Whiz- 
ard. After initializing FeynRules and loading a model 
definition, the code can be invoked in a similar fashion 
as the other interfaces described in Ref. [T3] in order to 
transform the Feynman rules obtained automatically from 
a Lagrangian into a set of model hies ready to use with 
Whizard. Both the legacy branch 1.9x as well as the cur- 
rent 2.x series of Whizard are actively supported by the 
code generator. Options that can be passed to the interface 
are introduced only as needed, while a detailed explana- 
tion of all options can be found in Appendix |XJ 



Supported fields and vertices: All helds currently sup- 
ported by FeynRules (scalars, Dirac and Majorana fer- 
mions, vectors and symmetric tensors) are handled by the 
interface. The set of accepted operators, the full list of 
which can be found in Tab. [1] is a subset of all the opera- 
tors supported by O'Mega. While still limited, this list is 
sufficient for a large number of BSM models. In addition, 
a future version of Whizard/ O'Mega will support the 
definition of completely general Lorentz structures in the 
model, allowing the interface to translate all interactions 
handled by FeynRules. 



Color: Color is treated in O'Mega in the color flow de- 
composition, with the flow structure being implicitly de- 
termined from the representations of the particles present 
at the vertex. Therefore, the interface has to strip the color 
structure from the vertices derived by FeynRules before 
writing them out to the model flies. While this process is 
straightforward for all color structures which correspond 
only to a single flow assignment, vertices with several pos- 
sible flow configurations must be treated with care in or- 
der to avoid mismatches between the flows assigned by 
O'Mega and those actually encoded in the couplings. To 
this end, the interface derives the color flow decomposi- 
tion from the color structure determined by FeynRules 
and rejects all vertices which would lead to a wrong flow 
assignment by O'Mega (these rejections are accompanied 
by warnings from the interface). 

At the moment, the SU(3)c representations supported 
by both Whizard and the interface are singlets (1), triplets 
(3), antitriplets (3) and octets (8). Tab. [2] shows all com- 
binations of these representations which can form singlets 
together with the support status of the respective color 
structures in Whizard and the interface. Although the 
supported color structures do not comprise all possible sin- 
glets, the list is sufficient for a large number of SM exten- 
sions. Furthermore, afuture revision of Whizard/ O'Mega 
will allow for explicit color flow assignments, thus remov- 
ing most of the current restrictions. 

In models generated for the 1.9x branch of Whizard, 
a hardcoded limit exists on the maximum number n c f of 
external color lines which is related to the number of ex- 
ternal octets n% and triplets / antitriplets 713 as 

n Gf = n 8 + — . 

This limit is set to n c { = 4 by default and can be changed 
via the option WOMaxNcf . No such limit applies to Whizard 
2.0.0 and higher. 



Running as- Starting with version 2.0, Whizard sup- 
ports a running strong coupling. While this is fully sup- 
ported by the interface, a choice has to be made which 
quantities are to be reevaluated when the strong coupling- 
is evolved. We choose to implement the running such that 
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Particle spins 


Supported Lorentz structures 


FFS 


All operators of dimension four are supported. 


FFV 


All operators of dimension four are supported. 


SSS 


All dimension three interactions are supported. 


svv 


Supported operators: 

dimension 3: O3 = V/'V^^ 

dimension 5- 0*. — 6 f<9 M W — &"V?\ (d„Vo„ — d„V->„\ 
Note that O5 generates the effective gluon-gluon-Higgs couplings obtained by 
integrating out heavy quarks. 


ssv 


(<^i<9 M <^>2 — 02<9''0i) Vfi type interactions are supported. 


ssvv 


All dimension four interactions are supported. 


ssss 


All dimension four interactions are supported. 


vvv 


All parity-conserving dimension four operators are supported, with the restriction that 
non-gauge interactions may be split into several vertices and can only be handled if all 
three fields are mutually different. 


WW 


All parity conserving dimension four operators are supported. 


TSS, TW, TFF 


The three point couplings in the Appendix of Ref. |45| are supported. 



of the particles. "S" stands for scalar, "F" for fermion (either Majorana or Dirac) and "V" for vector. 



SU(3)c representations 


Support status 


111, 331, 338, 
1111, 3311, 3381 


Fully supported by the interface 


888, 8881 


Supported only if at least two of the octets are identical particles. 


881, 8811 


Supported in output for Whizard 1.96 and higher (including the new 2.x branch) onlj0. 


3388 


Supported only if the octets are identical particles. 


8888 


The only supported flow structure is 
1 2 

^ ^ • r(l,2, 3, 4) + all acyclic permutations 
4 / ""*" X 3 

where f(l, 2, 3, 4) represents the Lorentz structure associated with the first flow. 


333, 333, 3331 
3331, 3333 


Unsupported. 



Table 2. All possible combinations of three or four SU(3)c representations supported by FeynRules which can be used to 
build singlets, together with the support status of the corresponding color structures in Whizard and the interface. 



by default aS, G (see Ref. [H] for the nomenclature regard- 
ing the QCD coupling) and any vertex factors depending 
on them are evolved. The list of internal parameters that 
are to be recalculated (together with the vertex factors 
depending on them) can be extended (beyond aS and G) 
by using the option WORunParameters when calling the 
interface. As older versions of Whizard do not evolve the 
strong coupling, the interface does not support running 
with Whizard 1.x. 



If the selected gauge is Feynman or , the interface 
can automatically assign the proper masses to the Gold- 
stone bosons. This behavior is requested by using the 
WOAutoGauge option. In the R^ gauges, the symbol repre- 
senting the gauge £ must be communicated to the interface 
by using the WOGaugeSymbol option (the symbol is auto- 
matically introduced into the list of external parameters 
if WOAutoGauge is selected at the same time). This feature 
can be used to automatically extend models implemented 
in Feynman gauge to the R^ gauges. 



Gauge choices: The interface supports the unitarity, 
Feynman and R^ gauges. The choice of gauge must be 
communicated to the interface via the option WOGauge. 
Note that massless gauge bosons are always treated in 
Feynman gauge. 



Since Whizard is a tree-level tool working with helic- 
ity amplitudes, the ghost sector is irrelevant for Whizard 
and can be left out. Ghosts and interactions involving 
ghosts are dropped by the interface. 
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3.3 Usage 

Installation and basic usage: From version 1.6.0 on- 
ward, the interface is included into the official FeynRules 
distribution. In addition, the latest version can be down- 
loaded from the Whizard homepage on HepForge. An 
installer is included which installs the interface into an 
existing FeynRules installation (this allows to use the 
program with FeynRules 1.4.x where it is not part of 
the package). 

Once installed, the interface can be called and used in 
the same way as the other interfaces described in Ref . [T^] . 
For example, once the FeynRules environment has been 
initialized and a model has been loaded, the command 

WriteWOOutput [L] 

will call FeynmanRules to extract the Feynman rules from 
the Lagrangian L, translate them together with the model 
data and finally write the files necessary for using the 
model with Whizard to an output directory (the name 
of which is inferred from the model name by default). Op- 
tions can be added for further control over the translation 
process (see Appendix . Instead of using a Lagrangian, 
it is also possible to call the interface on a vertex list. For 
example, the following command 

WriteWOOutput [Input -> list] 

will directly translate the vertex list list. Note that this 
vertex list must be given in flavor-expanded form in order 
for the interface to process it correctly. 

The interface also supports the WriteWOExtParams com- 
mand described in Ref. [12] for the other interfaces. Issuing 

WriteWOExtParams [filename] 

will write a list of all the external parameters to filename. 
The format of the output depends on the Whizard ver- 
sion it is targeted at: for 1.9x, it is a Fortran namelist 
suitable for use in whizard . in, while a Sindarin script 
is generated for 2.x. The target version is selected with 
the option WOWhizardVersion (which is the only option 
accepted by the command). 

During execution, the interface will print out a series 
of messages. It is highly advised to carefully read through 
this output as it not only summarizes the settings and the 
location of the output files, but also contains information 
on any skipped vertices or potential incompatibilities of 
the model with Whizard. 



Propagating the output to Whizard: After the in- 
terface has run successfully and written the model files 
to the output directory, the model must be imported into 
Whizard. This step depends on whether the code has 
been generated for the legacy 1.9x branch or for 2.x. 

In the case of version 1.9x, the model files must be 
patched into the Whizard tree before the model can be 
used. For this purpose, a script called inj ect . sh is created 
in the output directory. For most applications, this step is 
as easy as issuing (assuming that the output directory is 
the current working directory) 



. /inject. sh /path/to/whizard 

(replacing the path with the correct path to the Whizard 
installation). After reconfiguring and recompiling Whiz- 
ard, the model is available and can be used in a similar 
fashion as the built-in ones. The installer supports sev- 
eral options for more complex usage scenarios which are 
described in the INSTALL created in the output directory. 

For the new 2.x branch, the model files are compiled 
and can then be installed independently of Whizard. In 
the simplest scenario, assuming that the output directory 
is the current working directory and that the Whizard 
binaries can be found in the current ${PATH}, the instal- 
lation is performed by issuing 

./configure kk make clean kk make install 

This will compile the model and install it into the di- 
rectory ${H0ME}/ .whizard, making it fully available to 
Whizard without any further intervention. The build sys- 
tem can be adapted to more complicated cases through 
several options to configure which are listed in the INSTALL 
file created in the output directory. 

3.4 Validation of the interface 

The output of the interface has been extensively validated 
for both the 1.9x and 2.x branches. Specifically, the inte- 
grated cross sections for all possible 2 — » 2 processes in the 
FeynRules SM, the MSSM and the Three-Site Higgsless 
Model [13] have been compared between Whizard, Mad- 
Graph and CalcHep (for all three of which interfaces are 
included in FeynRules). In addition, the native versions 
of these models for the respective codes have also been 
included into the comparison^]. For those codes which sup- 
port it (CalcHep and Whizard), different gauges have 
been checked as well. In all comparisons, excellent agree- 
ment within the Monte Carlo errors was achieved. Several 
tables showcasing the results for selected processes in all 
three models can be found in Appendix [B] 



4 Usage examples 

In order to illustrate the usage of the interface as detailed 
above, we will now show several examples of generating 
output for Whizard 2. The examples are constructed to 
show the application of the different options of the inter- 
face and to serve as a starting point for the generation of 
Whizard versions of other FeynRules models. 



4.1 Standard Model 

To start off, we will create Whizard 2 versions of the 
Standard Model as implemented in FeynRules for dif- 
ferent gauge choices. 

4 A native version of the Three-Site Model is available only 
for CalcHEP and Whizard. 
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Unitarity Gauge In order to invoke FeynRules, we 
change to the corresponding directory and load the pro- 
gram in Mathematica via 

$FeynRulesPath = 

SetDirectory ["<path-to-FeynRules>"] ; 
<<FeynRules ' 

The model is loaded by 

LoadModel["Models/SM/SM.fr"] ; 
FeynmanGauge = False; 

Note that the second line is required to switch the Stan- 
dard Model to Unitarity gauge as opposed to Feynman 
gauge (which is the default). Generating a Whizard 2 
version of the model is now as easy as doing 

WriteWOOutput [LSM] ; 

After invokation, the interface first gives a short sum- 
mary of the setup 

Short model name is "f r_standard_model" 
Gauge: Unitarity 

Generating code for WHIZARD / O'Mega 

version 2.0.3 
Maximum number of couplings per FORTRAN 

module : 500 
Extensive lorentz structure checks disabled. 

Note that, as we have not changed any options, those set- 
tings represent the defaults. The output proceeds with 
the calculation of the Feynman rules from the Standard 
Model Lagrangian LSM. After the rules have been derived, 
the interface starts generating output and tries to match 
the vertices to their Whizard / O'Mega counterparts. 

10 of 75 vertices processed. . . 
20 of 75 vertices processed... 
30 of 75 vertices processed. . . 
40 of 75 vertices processed. . . 
50 of 75 vertices processed... 
60 of 75 vertices processed... 
70 of 75 vertices processed. . . 
processed a total of 75 vertices, kept 74 
of them and threw away 1 , 1 of which 
contained ghosts or goldstone bosons. 

The last line of the above output is particularily interest- 
ing, as it informs us that everything worked out correctly: 
the interface was able to match all vertices, and the only 
discarded vertex was the QCD ghost interaction. After the 
interface has finished running, the model files in the out- 
put directory are ready to use and can be compiled using 
the procedure described in the previous section. 



Feynman and R^ gauges As the Standard Model as im- 
plemented in FeynRules also supports Feynman gauge, 
we can use the program to generate a Feynman gauge 
version of the model. Loading FeynRules and the model 



proceeds as above, with the only difference being the change 
FeynmanGauge = True; 

In order to inform the interface about the modified gauge, 
we have to add the option WDGauge 

WriteWOOutput [LSM, WOGauge -> WOFeynman] ; 

The modified gauge is reflected in the output of the inter- 
face 

Short model name is "f r_standard_model" 
Gauge : Feynman 

Generating code for WHIZARD / O'Mega 

version 2.0.3 
Maximum number of couplings per FORTRAN 

module : 500 
Extensive lorentz structure checks disabled. 

The summary of the vertex identification now takes the 
following form 

processed a total of 163 vertices, kept 139 
of them and threw away 24, 24 of which 
contained ghosts . 

Again, this line tells us that there were no problems - 
the only discarded interactions involved the ghost sector 
which is irrelevant for Whizard. 

As Whizard does not depend on the ghost sector, the 
only difference between the different gauges from the per- 
spective of the interface are the gauge boson propagators 
and the Goldstone boson masses. Therefore, the interface 
can automatically convert a model in Feynman gauge to 
a model in gauge. To this end, the call to the interface 
must be changed to 

WriteWOOutput [LSM, WOGauge -> WORxi , 
WOAutoGauge -> True] ; 

The WOAutoGauge argument instructs the interface to au- 
tomatically 

1. Introduce a symbol for the gauge parameter £ into the 
list of external parameters 

2. Generate the Goldstone boson masses from those of the 
associated gauge bosons (ignoring the values provided 
by FeynRules) 

The modified setup is again reflected in the interface out- 
put 

Short model name is "f r_standard_model" 

Gauge : Rxi 

Gauge symbol: "Rxi" 

Generating code for WHIZARD / O'Mega 

version 2.0.3 
Maximum number of couplings per FORTRAN 

module : 500 
Extensive lorentz structure checks disabled. 

Note the default choice Rxi for the name of the £ parame- 
ter — this can be modified via the option WOGaugeParameter. 

While the WOAutoGauge feature allows to generate R^ 
gauged models from models implemented in Feynman gauge, 
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it is of course also possible to use models genuinely imple- 
mented in i?£ gauge by setting this parameter to False. 
Also, note that the choice of gauge only affects the propa- 
gators of massive fields. Massless gauge bosons are always 
treated in Feynman gauge. 



Compilation and usage In order to compile and use 
the freshly generated model files, change to the output 
directory which can be determined from the interface out- 
put (in this example, it is f r_standard_model-WO). As- 
suming that Whizard is available in the binary search 
path, compilation and installation proceeds as described 
in the previous Section by issuing (from the shell) 

./configure && make && make install 

The model is now ready and can be used similarly to 
the builtin Whizard models. For example, a minimal 
Whizard input file for calculating the e + e~ — ► W + W~ 
scattering cross section in the freshly generated model 
would look like 

model = f r_standard_model 

process test = "e+", "e-" -> "W+" , "W-" 

sqrts = 500 GeV 

integrate (test) 



4.2 Minimal Supersymmetric Standard Model 

In this Section, we illustrate the usage of the interface 
between FeynRules and Whizard in the context of the 
MSSM. For the conventions on the names of the particles 
and the relative signs entering the different terms of the 
Lagrangian, we refer to Refs. [T3lll5| . All the parameters 
of the model are then ordered in Les Houches blocks and 
counters following the SUSY Les Houches Accord (SLHA) 
[47,415]. Let us emphasize that implementing any model in 
FeynRules only requires a SLHA-like structure, whilst 
the choice of the SLHA conventions for the MSSM is only 
due to the author of the model implementation. 

The neutralino sector deserves special attention. After 
diagonalization of the mass matrix expresssed in terms of 
the gaugino and higgsino eigenstates, the resulting mass 
eigenvalues may be either negative or positive. In this case, 
two procedures can be followed. Either the masses are 
rendered positive and the associated mixing matrix gets 
purely imaginary entries or the masses are kept signed, 
the mixing matrix in this case being real. According to the 
SLHA agreement, the second option is adopted. For a spe- 
cific eigenvalue, the phase is absorbed into the definition 
of the relevant eigenvector, rendering the mass negative. 
However, Whizard does not officially support negative 
masses. For the case of the interface to FeynRules this 
means, that one must be careful using a SLHA file with 
explicit factors of the complex unity in the mixing ma- 
trix, and on the other hand, real and positive masses for 
the neutralinos. For the hard-coded SUSY models, this is 



completely handled internally, and the SUSY implemen- 
tations in WHIZARD have been extensively tested [23,25, 
26,33,34,35 . Especially Ref. [53] discusses the details of 
the neutralino (and chargino) mixing matrix. 

After having downloaded the model from the Feyn- 
Rules website, we store it in a new directory, labelled 
MSSM, of the model library of the local installation of Feyn- 
Rules. The model can then be loaded in Mathematica 
as 

$FeynRulesPath = 

SetDirectory ["<path-to-FeynRules>"] ; 
<<FeynRules ' 

LoadModel["Models/MSSM/MSSM.fr"] ; 
FeynmanGauge = False; 

As it can be seen from the last command line, we have 
adopted the unitarity gauge. Other choices, such as Feyn- 
man gauge, are also possible, and we refer to Section 14.11 
for using the interface with other gauge choice. 

The number of vertices associated to supersymmet- 
ric Lagrangians is in general very large (several thou- 
sands). For such models with many interactions, it is rec- 
ommended to first extract all the Feynman rules of the 
theory before calling the interface between FeynRules 
and Whizard. The reason is related to the efficiency 
of the interface which takes a lot of time in the extrac- 
tion of the interaction vertices. In the case one wishes to 
study the phenomenology of several benchmark scenar- 
ios, this procedure, which is illustrated below, allows to 
use the interface in the most optimal way. The Feynman 
rules are derived from the Lagrangian once and for all 
and then reused by the interface for each set of Whizard 
model files to be produced, considerably speeding up the 
generation of multiple model files issued from a single 
Lagrangian. In addition, the scalar potential of super- 
symmetric theories contains a large set of four scalar in- 
teractions, in general irrelevant for collider phenomenol- 
ogy. These vertices can be neglected with the help of the 
Exclude4Scalars->True option of both the commands 
FeynmanRules and WriteWODutput. The Feynman rules 
of the MSSM are then computed by issuing in the Math- 
ematica notebook, 

rules = FeynmanRules [lag, 

Exclude4Scalars->True , FlavorExpand->True] ; 

where lag is the variable containing the Lagrangian. 

By default, all the parameters of the model are set 
to the value of 1. A complete parameter file must there- 
fore be loaded. We choose the SPSlaE] benchmark point 
|49j . where we take care of setting all the masses of the 
neutralinos to positive values and adapt the associated 
mixing matrix accordingly. This parameter file, named 
spsla_wo.dat, can be downloaded from the FeynRules 
website, and loaded into FeynRules as 



Though we are aware that this specific point is excluded by 
LHC data, we use it here purely for demonstration purposes, 
as the validation of the MSSM within the interface has been 
done with it. 
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ReadLHAFile [Input -> "spsla_wo.dat"]; 

We note that this command does not reduce the size of 
the model output by removing vertices with vanishing cou- 
plings. However, if desired, this task could be done with 
the LoadRestriction command (see Ref. [SD] for details). 

The vertices are now ready to be exported to Whizard 
with the help of the command 

WriteWOOutput [Input -> rules] ; 

For compilation and usage of the produced model, we re- 
fer to Section B~T1 Finally, let us note that the numerical 
values of the parameters of the model can be modified 
directly at the level of the Whizard program, without 
having to generate a second time the Whizard model 
files from FeynRules. To this end, a Sindarin script is 
created by FeynRules with the help of the instruction 

WriteWOExtParams ["parameters . sin"] ; 

and can be further modified according to the needs of the 
user when employing the Whizard program. 



4.3 Three-Site Higgsless Model 

The Three-Site or Minimal Higgsless Model (MHM) gB] 
is a minimal deconstructed Higgsless model which con- 
tains only the first resonance in the tower of Kaluza-Klein 
modes of a Higgsless extra-dimensional model. It is a non- 
renormalizable, effective theory whose gauge group is an 
extension of the SM with an extra SU(2) gauge group. 
The breaking of the extended electroweak gauge symme- 
try is accomplished by a set of nonlinear sigma fields which 
represent the effects of physics at a higher scale and make 
the theory nonrenormalizable. The physical vector boson 
spectrum contains the usual photon, and Z bosons 
as well as a W 7± and Z' boson. Additionally, a new set 
of heavy fermions are introduced to accompany the new 
gauge group "site" which mix to form the physical eigen- 
states. This mixing is controlled by the small mixing pa- 
rameter €l which is adjusted to satisfy constraints from 
precision observables, such as the S parameter |51j . 

The MHM has been implemented into LanHEP [55], 
FeynRules [13] and independently into Whizard [53] . 
and the collider phenomenology has been studied by mak- 
ing use of these implementations [52,22,53 . Furthermore, 
the independent implementations in FeynRules and di- 
rectly into Whizard have been compared and found to 
agree (see Appendix [B]) . In this Section, we describe the 
use of the Whizard interface to generate Whizard model 
files for the MHM. This model has been implemented in 
Feynman gauge as well as unitarity gauge and contains the 
variable FeynmanGauge which can be set to True or False. 
When set to True, the option WDGauge-> WOFeynman must 
be used, as described in Appendix E] gauge can also 
be accomplished with this model by use of the options 
WDGauge -> WORxi and WOAutoGauge -> True. 



Since this model makes use of the nonlinear sigma field 

S = 1 + in - ^tt 2 H (1) 

many higher dimensional operators are included in the 
model which are not supported by Whizard. Although 
Whizard can reject these vertices and print a warning 
message to the user, it is preferable to remove the ver- 
tices with the option MaxCanonicalDimension->4 which 
is passed to FeynmanRules and restricts the Feynman rules 
to those of dimension four and sm 

Since the use of different gauges was already illustrated 
in Section 14.11 we will simply make use of this interface 
to output the model in Feynman gauge. We load Feyn- 
Rules with 

$FeynRulesPath = 

SetDirectory ["<path-to-FeynRules>"] ; 
<<FeynRules ' 

where <path-to-FeynRules> is the path to the Feyn- 
Rules files. The MHM is loaded by 

SetDirectory [ " <path-to-MHM> " ] ; 
LoadModel ["3-Site-particles . f r " , 

"3-Site-parameters . f r" , 

"3-Site-lagrangian. f r"] ; 
FeynmanGauge = True ; 

where <path-to-MHM> is the path to the directory where 
the MHM model files are stored and where the output of 
the Whizard interface will be written. The Whizard in- 
terface is initiated with the command 

WriteWOOutput [LGauge , LGold, LGhost, LFermion, 
LGoldLeptons , LGoldQuarks, 
MaxCanonicalDimension->4 , 

WOGauge->WOFeynman, WOModelName->"f r_mhm"] ; 

where we have also made use of the option WOModelName 
to change the name of the model as seen by Whizard (see 
Appendix [SJ . As in the SM (see Section l4"7Tj) . the interface 
begins by writing a short informational message: 

Short model name is "fr_mhm" 
Gauge : Feynman 

Generating code for WHIZARD / O'Mega 

version 2.0.3 
Automagically assigning Goldstone 

boson masses . . . 
Maximum number of couplings per FORTRAN 

module : 500 
Extensive lorentz structure checks disabled. 

After calculating the Feynman rules and processing the 
vertices, the interface gives the summary: 



6 We note, at this stage, that MaxCanonicalDimension is an 
option of the FeynmanRules function rather than of the in- 
terface, itself. In fact, the interface accepts all the options of 
FeynmanRules and simply passes them on to the latter. 
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processed a total of 922 vertices, kept 633 
of them and threw away 289, 289 of which 
contained ghosts. 

showing that no vertices were missed. The files are stored 
in the directory f r_mhm and are ready to be installed and 
used with Whizard. 



5 Conclusions 

In this paper, we presented the interface between Feyn- 
Rules, a program to automatically generate Feynman 
rules from a generic input Lagrangian, to the multi-particle 
event generator Whizard. The interface can be used both 
for the legacy version Whizard l.xx as well as the actual 
release line Whizard 2.x and is included in FeynRules 
from version 1.6 on. There has been a quite exhaustive 
validation of the interface both with FeynRules mod- 
els used with other Monte Carlo generators as well as 
with hard-coded implementations of models. We tried to 
keep the description of the interface rather self-contained, 
but recommend the user to consult both Whizard and 
FeynRules manuals when performing a calculation with 
the interface. The paper gives a complete overview of the 
supported particle types as well as the Lorentz and color 
structures of interaction vertices. We listed the options 
and flags which steer the compilation of physics models 
and the transferring of the Feynman rules into the gener- 
ator. In order to give examples of the workflow, we applied 
it to the three models used in our validation procedure. 
From the steps given in the paper to go from the La- 
grangian formulation of that model through FeynRules 
to a simulation within Whizard, the user can easily ap- 
ply the interface to his or her own preferred model. The 
presented FeynRules/ Whizard interface hopefully fa- 
cilitates phenomcnological studies at present and future 
collider experiments. A prime example of using the com- 
bination between FeynRules and an event generator like 
Whizard is to investigate models whose particle content 
depend on the point in parameter space, which is however 
a separate project and will be presented elsewhere [54] . 
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WOWhizardVersion 


Whizard versions supported 


"1.92" 


1.92 


"1.93" 


1.93 - 1.95 


"1.96" 


1.96+ 


"2.0" 


2.0.0 - 2.0.2 


"2.0.3" (default) 


2.0.3+ 



Table 3. Currently available choices for the 
WOWhizardVersion option, together with the respective 
Whizard versions supported by them. 



A Options of the Whizard interface 

In the following we present a comprehensive list of all 
the options accepted by WriteWOOutput. Additionally, we 
note that all options of FeynmanRules are accepted by 
WriteWOOutput, which passes them on to FeynmanRules. 

Input 

An optional vertex list to use instead of a Lagrangian 
(which can then be omitted). 

WOWhizardVersion 

Select the Whizard version for which code is to be 
generated. The currently available choices are summa- 
rized in Tab. [3] This list will expand as the program 
evolves. To get a summary of all choices available in 
a particular version of the interface, use the command 
?W0WhizardVersion. 

WOModelName 

The name under which the model will be known to 
Whizard. Please note that models for versions 1.9x 
must start with "fr_" if they are to be picked up 
by Whizard automatically. The default is determined 
from the FeynRules model name. 
Output 

The name of the output directory. The default is de- 
termined from the FeynRules model name. 
W0 Gauge 

Gauge choice (c/. Section l3~2l) . Possible values are: 
WOFeynman, WORxi and WOUnitarity (default). 

WOGaugeParameter 

The external or internal parameter representing the 
gauge £ in the gauges (c/. Section I3.2[) . Default: 
Rxi 

WOAutoGauge 

Automatically assign the Goldstone boson masses in 
the Feynman and gauges and automatically append 
the symbol for £ to the parameter list in the R^ gauges. 
Default: False 
WOMaxNcf 

Maximum number of color flows provided for in the 
output. Essentially, this is n$ — \n-s where n$ is the 
maximum number of external color octets and 713 is the 
maximum number of external triplets and antitriplets. 
Note that this only affects Whizard 1.9x; there is no 
limit for 2.0.0+ . Default: 4 
WORunParameters 

The list of all internal parameters which will be recal- 
culated if as is evolved (see above). This only affects 
Whizard 2.0.0+. Default: {aS, G} 
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WDFast 

If the interface drops vertices which are supported, this 
option can be set to False to enable some more time 
consuming checks which might aid the identification. 
Default: True 

WDMaxCouplingsPerFile 

The maximum number of couplings that are written 
to a single Fortran file. If compilation takes too long 
or fails, this can be lowered. Default: 500 

WDVerbose 

Enable verbose output and in particular more exten- 
sive information on any skipped vertices. Default: False 

B Validation Tables 

In this Appendix we summarize the validation of the in- 
terface from FeynRules to Whizard which was per- 
formed by calculating the integrated cross sections related 
to all possible two-to-two partonic processes in the SM, 
the MSSM and the Minimal Higgsless Model. We evalu- 
ated this quantity both with the FEYNRuLES-generated 
model files and the built-in ("stock") implementations of 
the various matrix-element generators, if existing, in the 
framework of CalcHep, MadGraph/MadEvent and 
Whizard/O'mega. Where supported by both model and 
generator, we also checked both Feynman and unitarity 
gauge. We ran the codes such that the Monte Carlo uncer- 
tainty was pushed below one percent, and excellent agree- 
ment at this level was achieved in all cases. In Tab. Eland [7] 
we present a selection of the results we obtained. The full 
list of results is available from the FeynRules website. 

In order to avoid collinear divergences, we imposed a 
cut on the px of each final state particle, 



For the MSSM, the parameter set used corresponds to the 
benchmark point SPSla (with the strong coupling set 
to a s (M z ) = 0.1172), whereas for the SM and MHM the 
values of the input parameters are summarized in Tab. U 
and [3] The strong coupling was evolved to y/s for each 
process. All particle widths were set to zero. 



Parameter 


Symbol 


Value 


Inverse of the 






electromagnetic coupling 


„ —1 f Ti/r \ 
a Ew( M Z) 


127.9 


Strong coupling 


a s {Mz) 


0.1172 


Z mass 


Mz 


91.1875 Gev 


W mass 


IVlw 


on one r*r.\T 


W' mass 


M W i 


500 GeV 


Heavy fermion mass 


M F 


4TeV 


H lepton mass 


m M 


0.1057 GeV 


t lepton mass 


m T 


1.777 GeV 


s quark mass 


m s 


0.104 GeV 


c quark mass 


m c 


1.27 GeV 


b quark mass 


m b 


4.2 GeV 


t quark mass 


m t 


171.2 GeV 



Table 5. Input parameters for the MHM. 



Parameter 


Symbol 


Value 


Inverse of the 






electromagnetic coupling 




127.9 


Strong coupling 


a s {M z ) 


0.1172 


Fermi constant 


G F 


1.16639e-5 GeV -2 


Z pole mass 


M z 


91.188 GeV 


c quark mass 


m c 





b quark mass 


m b 


4.7 GeV 


t quark mass 


mt 


174.3 GeV 


t lepton mass 


m T 


1.777 GeV 


Higgs mass 


m H 


120 GeV 


Cabibbo angle 


6c 


0.227736 



Table 4. Input parameters for the SM. 
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0.250 


0.250 


0.250 


w + w~ 


-4 


H H 


1599 


0.151 


0.151 


0.151 


0.151 


0.151 


0.151 


0.151 


0.151 


0.151 


0.151 


t t 


— y 


H H 


2354 


0.0593 


0.0593 


0.0593 


0.0593 


0.0593 


0.0593 


0.0593 


0.0593 


0.0593 


0.0593 




— >■ 


H H 


974 


2.96e-06 


2.96e-06 


2.96e-06 


2.96e-06 


2.96e-06 


2.96e-06 


2.96e-06 


2.96c-06 


2.96e-06 


2.96e-06 


b b 


— ► 


H H 


998 


6.33c-06 


6.33c-06 


6.33c-06 


6.33c-06 


6.33o-06 


6.33o-06 


6.33c-06 


6.33c-06 


6.33c-06 


6.33c-06 



Table 6. Cross sections for a selection of two-to-two processes in the SM. The results for CalcHep, MadGraph /MadEvent, 
Whizard 1.x and Whizard 2.0 are labeled by CH, MG, WOl and W02, respectively. 'F' and 'LP refer to the choice of the 
gauge in which the calculation was performed, Feyman or unitarity. The pr cuts were fixed according to Eq. and the unit 
for all cross sections is pb. 
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FeynRules 








Stock 




Process 


V a 


MG 


CH 


WO 2 


WOl 


WOl 


CH 


W02 


W02 








fGcVl 


U 


U 


u 


F 


U 


F 


F 


u 


9 9 


"* 


9 9 


200 


1.88c04 


1.88c04 


1.88c04 


1.88c04 


1.88c04 


1.88c04 


1.88e04 


1.88e04 


Z W~i~ 


— y 


Z 1 W' 


4693 


103 


103 


103 


103 


103 


103 


103 


103 


7 7 


— y 


W~t \\~ ~ 


643 


16. 1 


16. 1 


16.1 


16. 1 


16.1 


16. 1 


16.1 


16. 1 


z 


— ¥ 


W' z' 


3015 


0.755 


0.755 


0.755 


0.755 


0.755 


0.755 


0.755 


0.756 


7 Z 


— > 


w~^~ w~ 


8730 


0.0412 


0.0412 


0.0411 


0.0412 


0.0412 


0.0412 


0.0412 


0.0412 


7 Z 


— > 


W w~ 


2686 


1.30 


1.30 


1.30 


1.30 


1.30 


1.30 


1.30 


1.30 




— y 




6322 


3.32 


3.31 


3.31 


3.31 


3.31 


3.31 


3.31 


3.32 


z z' 


— y 


w~^~ w~ 


6371 


3.86 


3.86 


3.86 


3.86 


3.86 


3.86 


3.86 


3.86 








2965 


2.41 


2.42 


2.42 


2.42 


2.42 


2.42 


2.42 


2.42 


7 W T 




z' 


8656 


0.0349 


0.0349 


n n o A n 

0.0349 


0.0349 


0.0349 


o.034y 


n n o A n 

0.0349 


0.034y 






s s' 


18765 


33.6 


33.6 


33.6 


33.6 


33.6 


33.6 


33.6 


33.6 


z' z' 


— y 


b b 


8094 


0.0958 


0.0959 


0.0959 


0.0959 


0.0959 


0.0959 


0.0959 


0.0959 


z' z' 


— y 


— ' +' 


36900 


1.54 


1.54 


1.54 


1.54 


1.54 


1.54 


1.54 


1.54 




— y 


c' s' 


34886 


0.00220 


0.00220 


0.00220 


0.00220 


0.00220 


0.00220 


0.00220 


0.00220 




— y 


T T~^~ 


35208 


1.45 


1.45 


1.45 


1.45 


1.45 


1.45 


1.45 


1.45 


\Y+ W~ 


— y 


U U 


18765 


34.0 


34.0 


34.0 


34.0 


34.0 


34.0 


34.0 


34.0 


z z 


— y 


u u 


17173 


30.3 


30.3 


30.3 


30.3 


30.3 


30.3 


30.3 


30.3 


w~ 


—y- 


b b' 


20460 


38.5 


38.5 


38.5 


38.5 


38.5 


38.5 


38.5 


38.4 


w~^~ w~ 


— y 


v e ' v„ 


36886 


1.64 


1.64 


1.64 


1.64 


1.64 


1.64 


1.64 


1.64 


7 7 




(j, fi+ 


32886 


0.000463 


0.000463 


0.000463 


0.000463 


0.000463 


0.000463 


0.000463 


0.000463 


V-T u 




u T u' 


32886 


0.0966 


0.0966 


0.0967 


0.0966 


0.0966 


0.0966 


0.0966 


0.0966 


I/ e u 


— > 


u e ' u 


32886 


u.uyoo 


u.uyoo 


u.uyoo 


u.uyou 


u.uyoo 


u.uyou 


0.0967 


0.0967 


c c 


— y 




32896 


5.38c-06 


5.38e-06 


5.38e-06 


5.38e-06 


5.38e-06 


5.38e-06 


5.38e-06 


5.38o-06 


b u' 


—y- 


t' d' 


49606 


0.0337 


0.0337 


0.0337 


0.0337 


0.0337 


0.0337 


0.0337 


0.0337 


r~ d 


—y- 


V T ' u 


32893 


0.386 


0.386 


0.386 


0.386 


0.386 


0.386 


0.386 


0.386 


e~ b 


— y 


e~' b' 


32903 


0.0965 


0.0965 


0.0965 


0.0965 


0.0965 


0.0965 


0.0966 


0.0965 


t u' r 


—y- 


T+' b' 


50014 


0.0430 


0.0429 


0.0429 


0.0429 


0.0429 


0.0429 


0.0429 


0.0429 


c s 


—y- 


u e ' 


32892 


1.82c-05 


1.82e-05 


1.82e-05 


1.82e-05 


1.826-05 


1.82e-05 


1.82e-05 


1.82c-05 




—y- 


t t' 


17388 


7.23c-06 


7.24e-06 


7.24e-06 


7.24e-06 


7.24e-06 


7.24e-06 


7.24e-06 


7.24c-06 


T~ T + 




s s 


16458 


9.33c-06 


9.35c-06 


9.35c-06 


9.35c-06 


9.35c-06 


9.35c-06 


9.35c-06 


9.35o-06 



Table 7. Cross sections for a selection of two-to-two processes in the MHM. The results for CalcHep, MadGraph /MadEvent, 
Whizard 1.x and Whizard 2.0 are denoted by CH, MG, WOl and W02, respectively. *F' and 'U' refer to the choice of the 
gauge in which the calculation was performed, Feyman or unitarity. The pr cuts were fixed according to Eq. and the unit 
for all cross sections is pb. A trailing prime "/" at the end of a particle identifier denotes its heavy partner (c/. Ref. [46]). 
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Process 



z z 

w + w + 
7 z 
a a 
z z 



7° v+ 
X4 Xi 



c c 
b b 



c c 
s s 



T U 
"fl V t 

7 z 
7 z 



w+ w~ 

w+ W+ 
W+ w~ 
a a 
z z 

w+ w~ 



a a 
7 z 
7 w+ 
z z 
a a 
7 7 

w+ w 
W+ w~ 



a a 



l 6 u 5 
Si l>2 

r i+ 

v 3 l~ 
U4 H~ 
"1 



[GeV] 



1368 

1277 

1003 

800 

1459 

1278 



4862 

758 

2573 

730 

1400 

800 

678 

1412 



1114 

2897 
5286 
2179 
4862 
1449 
3694 
4862 



3079 
1480 
1730 
3349 
4854 
1317 
4263 
3956 



FeynRules 
CH MG WOl W02 

u u u u 



26.2 

25.6 

19.3 

839 

0.223 

4.95 



26.3 

25.7 

19.3 

840 

0.223 

4.96 



26.2 
25.6 
19.3 

839 

0.223 

4.95 



26.2 
25.7 
19.3 

839 

0.223 

4.95 



0.977 

0.327 

0.0219 

0.0838 

6.00 

0.0272 

1.72 

0.00586 



0.977 

0.327 

0.0219 

0.0839 

6.00 

0.0272 

1.72 

0.00585 



0.977 

0.327 

0.0219 

0.0838 

6.00 

0.0272 

1.72 

0.00586 



0.977 

0.327 

0.0219 

0.0838 

5.99 

0.0272 

1.72 

0.00586 



0.000873 

0.00237 

0.0151 

0.000134 

0.144 

0.0291 

0.00759 

0.144 



0.000874 

0.00237 

0.0151 

0.000134 

0.144 

0.0291 

0.00759 

0.144 



0.000873 

0.00237 

0.0151 

0.000134 

0.144 

0.0291 

0.00758 

0.144 



0.000873 

0.00237 

0.0151 

0.000134 

0.144 

0.0291 

0.00759 

0.144 



0.000510 

0.00902 

0.00278 

9.85e-06 

0.00268 

0.0274 

9.84e-06 

0.000449 



0.000510 

0.00903 

0.00278 

9.85c-06 

0.00268 

0.0274 

9.84c-06 

0.000449 



0.000510 

0.00902 

0.00278 

9.85c-06 

0.00268 

0.0274 

9.84e-06 

0.000449 



0.000510 

0.00902 

0.00278 

9.85o-06 

0.00268 

0.0274 

9.84c-06 

0.000449 



Stock 
MG W02 

u u 



26.3 

25.7 

19.3 

839 

0.223 

4.96 



0.977 

0.327 

0.0219 

0.0839 

6.00 

0.0272 

1.72 

0.00587 



0.000874 

0.00236 

0.0151 

0.000134 

0.144 

0.0291 

0.00759 

0.144 



0.000510 

0.00903 

0.00278 

9.85c-06 

0.00268 

0.0274 

9.84e-06 

0.000449 



26.2 
25.6 
19.3 

839 

0.223 

4.95 



0.977 

0.327 

0.0219 

0.0839 

5.99 

0.0272 

1.72 

0.00586 



0.000873 

0.00237 

0.0151 

0.000134 

0.144 

0.0291 

0.00759 

0.144 



0.000510 

0.00902 

0.00278 

9.85o-06 

0.00268 

0.0274 

9.83e-06 

0.000449 



Table 8. Cross sections for a selection of two-to-two processes in the MSSM. The results for CalcHep, MadGraph /MadE- 
vent, Whizard 1.x and Whizard 2.0 are denoted by CH, MG, WOl and W02, respectively. 'F' and 'U' refer to the choice of 
the gauge in which the calculation was performed, Feyman or unitarity. The pr cuts were fixed according to Eq. l[2]l. and the 
unit for all cross sections is pb. The sfermions are labeled according to the SLHA2 conventions (cf. Ref. |48| ) . 
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